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Photonic quantum information processing schemes, such as linear optics quantum computing, 
and other experiments relying on single-photon interference, inherently require complete photon 
indistinguishability to enable the desired photonic interactions to take place. Mode-mismatch is the 
dominant cause of photon distinguishability in optical circuits. Here we study the effects of photon 
wave-packet shape on tolerance against the effects of mode-mismatch in linear optical circuits, and 
show that Gaussian distributed photons with large bandwidth are optimal. The result is general and 
holds for arbitrary linear optical circuits, including ones which allow for post-selection and classical 
feed-forward. 
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I. INTRODUCTION 



The interference of single photons plays a central role 
in linear optics quantum computing [l| and many other 
quantum optics experiments. One of the major obstacles 
facing experimentalists is mode-mismatch, whereby pho- 
ton indistinguishability is compromised within a circuit, 
resulting in loss of quantum gate fidelity. The effects 
of input distinguishability, and, more generally, mode- 
mismatch, in linear optical circuits has been studied in a 
number of situations 0, 0, • 

The tolerance of optical circuits against the effects 
of mode-mismatch is highly dependent upon the shape 
of the interacting photons' wave-packets. In this paper 
we study this effect and derive conditions which maxi- 
mize the tolerance of optical circuits against the effects 
of mode-mismatch. We begin by considering the most 
trivial case of two photons interacting on a beamsplit- 
ter. We then generalize our findings and show that for 
arbitrary linear optical circuits, including ones which in- 
corporate post-selection and classical feed-forward, toler- 
ance against the effects of mode-mismatch is maximized 
when utilizing Gaussian shaped photons, which are as 
broad as possible in the degrees of freedom in which 
mode-mismatch is introduced. 

Understanding the influence wave-packet shape has on 
the effects of mode-mismatch is important from a practi- 
cal perspective, where experimentalists must choose the 
most appropriate photon engineering techniques. We 
provide a discussion of such techniques, in the context 
of our findings, in Section llVl 



II. PROOF THAT GAUSSIAN IS OPTIMAL 
FOR TWO PHOTONS INTERACTING ON A 
BEAMSPLITTER 

We begin by considering the most basic linear op- 
tics network: two photon interference on a beamsplit- 
ter. A beamsplitter with reflectivity r/ is described by 
the Heisenberg equations of motion 



a ou t = Vrj& + a/1 — 77b 
bout = V?/t) - y/i - rjk 



(1) 



where we assume the phase-asymmetric beamsplitter 
convention, a and b are the usual photon annihilation 
operators for the two spatial modes. 

If we consider an r\ = 0.5 beamsplitter with a single 
photon incident at each input (i.e. |V>in) = |l)a|l)b)> the 
output state it given by 



IV>o 



^(|2)a|0) b 



|0)a|2) t 



(2) 



Thus we see complete suppression of single photon terms 
as a result of quantum interference. 

Next we consider the non-ideal case, where mode- 
mismatch is present. We model mode-mismatch in the 
same manner described by First we associate a 

wave-function ip{ x )i with input photons, where x is 
some photonic degree of freedom. Note that modeling 
mode- mismatch in a single degree of freedom is sufficient 
to characterize arbitrary mode-matching effects, due to 
the inherent indistinguishability of the effects of mode- 
mismatch in different degrees of freedom. Thus, photons 
are represented as 



I '/'photon) 



ip(x)eJ(x) dx\0) 



(3) 



where a}(x) is the photonic creation operator at the in- 
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this form is represented by displacing the photons' wave- 
functions. This has the effect of transforming the wave- 
function according to ip(x) — > ip(x — t), where r is the 
displacement parameter. 

If we input a photon into each input of an rj = 0.5 
beamsplitter, where there is a relative displacement be- 
tween the photons, we no longer observe complete sup- 
pression of the single photon terms as before. Instead 
the probability of measuring a coincidence between the 
outputs (i.e. a single photon at each output) is given by 



M) = \-\ 



ip(x)*ip(x — t) dt 



(4) 



For a derivation see Q. This expression has the prop- 
erty that < (N a N h ) < 0.5, with (N a N h ) = if and 
only if r = 0. This is completely equivalent to the ef- 
fect observed by Hong, Ou and Mandel (HOM) 0. The 
vanishing of (iV a iVb) at r = is widely referred to as the 
'HOM-dip'. 

The behavior of (N a Nb) for non-zero r is highly de- 
pendent on the form of i(>(x). We are therefore moti- 
vated to ask what form of ifi(x) minimizes the effect of r 
on (iVaiVb), i.e. which maximizes the system's tolerance 
against the effects of mode-mismatch. 

We assume that mode-mismatch is some unknown 
small deviation from the ideal case. This assumption is 
justified, since if the displacement parameter were known 
and large it could be corrected for. 

We ask what form of ip(x) minimizes the curvature of 
the function (iV a iVb) against r. Thus we aim to minimize 
the function 



S = 



d 2 



dr 2 



(N a N b ) 
\ — l ip(x)*ip (x)dx 



T = 

oo 



subject to the normalization constraint 
\i(;(x)\ 2 dx = 1 



(5) 



(6) 



This optimization has a trivial solution. Namely, if 
we choose any form for i/)(x) and let its width approach 
infinity, the function S will exhibit no dependence on r. 
Thus we can immediately establish the following criteria 
for optimal photon engineering: photon wave-functions 
should be as broad as possible in the degree of freedom 
in which mode-mismatch is introduced. For example, 
in the presence of temporal mode-mismatch we should 
make photons as temporally broad as possible. We note 
that this leads to a trade-off between temporal mode- 
mismatch and clock-speed: the longer photons are, the 
slower a circuit can be operated. Similarly, there are 
limitations in how broad photons can be in other degrees 
of freedom, for example spatially. This result is very 
intuitive, since we expect that the closer a wave-function 
is to being invariant under translations (i.e. broader), 



the more stable circuit operation will be against such 
translations. 

Next we impose the additional constraint that tp(x) 
must have fixed variance, 



l \ip(x)\ 2 dx = Ax 



(7) 



We have assumed the mean, x\ip(x)\ 2 , vanishes. 
This is allowed because our analysis is invariant under 
global translations in the x coordinate. Fixing the vari- 
ance allows us to avoid the trivial solution, and also 
gives us a means by which to compare different func- 
tions. Thus, the optimizing function ip( x ) corresponds 
to the function, which, for a given bandwidth, maximizes 
tolerance against the effects of mode-mismatch. 

Figure ^ illustrates the behavior of (N a Nb) against 
r for Lorentzian, double-sided Lorentzian and Gaussian 
wave- functions 1 , where Ax = 1. It is immediately obvi- 
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FIG. 1: Hong-Ou-Mandel dip for Gaussian (solid), Lorentzian 
(broadly dashed) and double-sided Lorentzian (finely dashed) 
wave-functions, where Ax = 1. 

ous from Figure^that wave-packet shape plays a critical 
role in tolerance against the effects of mode-mismatch. 
Most notably, for Lorentzian distributed photons, which 
are characterized by a discontinuity, the coincidence rate 
increases far more rapidly than either the Gaussian or 
double-sided Lorentzian cases, both of which are smooth 
functions. 

We can reformulate the optimization problem in a 
more familiar quantum mechanical description as 



5« W 2 |V> 
subject to the constraints 

W#) = i 



(8) 



(9) 
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1 Lorentzian: 4>( x ) = 
i>( x ) = \J 3§T i+§£T. Gaussian: tp{x) = f 



doublc-sidcd Lorentzian: 
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FIG. 2: An arbitrary linear optics network consisting of n 
single photon inputs, all characterized by wave-function ip(x), 
and m vacuum inputs. The detectors following some of the 
outputs facilitate post-selection. 



FIG. 3: Tk,i represents the displacement introduced into the 
photon wave-packet when it travels from the fcth input to the 
lih output. 



and 



= Aa 



(10) 



where x and p can be considered the usual position and 
momentum operators (but could equally correspond to 
any Fourier pair). The position-momentum Heisenberg 
uncertainty relation is given by 



((Axf)((Apf)> 1 

Upon applying the constraint ((Ax) 2 ) = 
reduces to 

((Ap) 2 )> 



(11) 



Ax, this 



4Aa 



(12) 



It is known that we obtain equality for Gaussian ij){x), 
i.e. a Gaussian state is a minimum uncertainty state [jj. 
However ((Ap) 2 ) = (p 2 ) oc S. Thus, for Gaussian ip(x), 
S is minimized, as required. 



III. GENERAL PROOF FOR ARBITRARY 
LINEAR OPTICS NETWORK 

We now consider an arbitrary linear optics network of 
the form shown in Figure [2] The input state, according 
to this model, is of the form 



\4>in) = / ip(x)a[(x) dx / ip(x)&2(x) dx . . . 

J — oo J — oo 

/oo 
V<z)4(z)dx|0) 
-oo 

: J] / ^(aOa}(*)da:|0) (13) 

where a|(x) are the photonic creation operators of the 
ith input. Circuits where inputs contain higher photon 
number terms are allowed for in this model by assum- 
ing a suitable beamsplitter network to be inside the box, 
followed by appropriate post-selection. 



If we allow U to be an arbitrary beamsplitter network 
acting on the input state, the output state will be of the 
form 



E 



n 



^{x)k\ 3 {x) dx\0) 
(14) 

which is simply a sum-of-paths of all possible routes the 
input photons (J) could take to reach all possible output 
configurations . . . , ij). The complex coefficients A 

are the amplitudes of particular paths through the cir- 
cuit. These parameters are a function of the circuit and 
completely characterize the output state. The A are ob- 
tained by tracing along paths from inputs to outputs. 
Upon reflection from or transmission through a beam- 
splitter, the respective parameter gains a factor equal to 
the beamsplitter's reflectivity or transmissivity respec- 
tively. Upon a phase change the parameter gains a com- 
plex rotation factor. 

Post-selection is accommodated for through suitable 
adjustment of the A parameters and discarding the de- 
grees of freedom associated with the measured modes. 
Classical feed-forward is accommodated for by recogniz- 
ing that a circuit with feed-forward can be broken down 
into multiple blocks of the form shown in Figure|21 where 
the A parameters in later blocks are determined by mea- 
surement outcomes from earlier blocks. 

We model mode-mismatch by introducing displace- 
ments into photon wave-packets as they travel between 
different inputs and outputs. We introduce the parame- 
ters Tk.i, which represent the cumulative displacement in- 
troduced between the fcth input and Zth output, shown in 
Figure[3] Thus, when a photon travels from the fcth input 
to the Zth output, its wave- function undergoes the trans- 
formation ip( x ) ~ > VK 2 - — T k,i)- We have assumed that 
input photons are all initially indistinguishable. There is 
no loss of generality in making this assumption since we 
can take the displacement parameters describing input 
distinguishability to be implicitly incorporated into the 
box. 

The general form of the output state in the presence 
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of mode-mismatch is given by 

IVWt) = E X ii,i2,-i n 

»1,«2 *i» = l 



We define the fidelity as the overlap between the ideal 
and non- ideal output states, 



71 />oo 

n / ^-^Kw^i ) ( i5 ) 



F = |<^out|^out>| S 



°l( E 'V ,, II / V^)*a, 3 (x)dx)( ^ A illi2 ,... in J] / #n-r^a].(z)dz)|0; 

n+rn n+m n ! 

E E V,. ; II / VW(z-r i3 . i4 )c 



(16) 



where the parameters /3 have been introduced to allow 
for the different combinations in which terms from the 
left-hand product can act on terms from the right-hand 
product. 

We apply the same criteria as for the HOM case, and 
attempt to find the form of ip(x) such that the curvature 
of the fidelity function is minimized, subject to the same 
normalization and variance constraints as before. Thus, 

for any given {m,n}. 



we wish to minimize 
The result is of the form 



8 2 F 



dr 2 



xjj(x)*ip (x) dx 



(17) 



which is the same minimization as previously. The pro- 
portionality factor is circuit dependent. Therefore, for an 
arbitrary linear optics network, the fidelity of the output 
state will be most resilient against the effects of mode- 
mismatch when input photons have a Gaussian profile 
in the degree of freedom in which the mode-mismatch is 
introduced. 



IV. DISCUSSION 



We now discuss the potential of various photon sources 
for quantum information processing in light of our re- 
sults. For the most part we assume cavity based sources 
producing beams with Gaussian spatial profiles and so 
concentrate on their temporal profiles, although we note 
that this is not necessarily the case for all experimental 
examples cited. 



A. Intra-cavity spontaneous photon emission 

When a 'fast' single photon emitter is placed in a 'slow' 
optical cavity we observe photons with an approximately 
Lorentzian frequency wave-function, 



K 



1 



7T K + lUJ 



(18) 



where u> is frequency and n is the cavity bandwidth. Here 
we are assuming that 7>k, where 7 is the spontaneous 
emission lifetime of the emitter. Examples of such sources 
include quantum dot 0,0 and fluorescence based 
sources. Typically such sources exhibit some inhomogc- 
nous broadening of the photon emission process, referred 
to as time-jitter. This results in a mixing effect whereby 
photons are characterized by a mixture of temporally 
displaced Lorentzian wave-functions. To minimize this 
effect the time-uncertainty of photon emission must be 
kept small compared to the decay time of the cavity. In 
the presence of time-jitter the state can be expressed in 
the form 



/(r) 



nOO pOO 
' — OO J —OO 



^ T il>{u)il}*{J)\uj)(J\dLodJ dr 

(19) 

where f(r) is determined by the time-jitter and charac- 
terizes the mixture. 

The effect of time-jitter on the Knill controlled-sign 
gate 0| has been examined by Kiraz et al. 0. It was 
found that to achieve gate fidelity of 99%, time-jitter 
must be kept below 0.3% of the inverse bandwidth. This 
is in stark contrast to an analysis of the simplified KLM 
CNOT gate [l3j indicating that, for Gaussian photons, 
gate fidelity of 99% requires temporal synchronization to 
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within 10% of the inverse bandwidth Jst]. This relative 
intolerance against time-jitter can be attributed to the 
temporal discontinuity inherent in the Lorenztian func- 
tion, and is the same reason we observe rapid falloff in 
the HOM dip for the Lorentzian case in Figure ^ 

We expect that photon sources based on spontaneous 
emission, and other sources which produce Lorentzian 
photons, will not be well suited to quantum information 
processing applications, unless suitable filtering or other 
shaping techniques are first applied. Most importantly, 
such techniques would have to eliminate the temporal dis- 
continuity inherent in the Lorentzian, which is the pri- 
mary culprit in loss of gate fidelity in the presence of 
temporal mode-mismatch. In their favor, however, such 
sources allow production of single photons on demand, 
which is very desirable for quantum information process- 
ing applications. 




FIG. 4: Conditional production of filtered photons through 
parametric down- conversion. The down-converter produces 
an entangled pair of photons into two distinct spatial modes. 
One of the modes is filtered and a photo-detector conditions 
upon the detection of a single photon following the filter. 
When post-selection succeeds, a photon will be present in 
the other path, whose wave-function will reflect that of the 
detected filtered photon. 



B. Non-degenerate parametric down-conversion 

Non-degenerate parametric down-conversion is widely 
used in quantum optics experiments for the production 
of heralded single photons. The down-conversion process 
probabilistically produces entangled photon pairs in dis- 
tinct spatial modes. The output from a down-converter 
in an optical cavity can be expressed in the form 

2 X k 



l^out) = |0) a |0) b + / „, 2 (|(U-„)a|UU,)b 

(20) 



+|l w 0_ w ) a |0 w l_ a) ) b )do; 



where a and b denote the spatial modes, \ is related to 
the conversion efficiency of the down-conversion process, 
and k to the bandwidth. It has been assumed that con- 
version efficiency is very weak, such that \ <C k, which 
justifies neglecting higher-order photon number terms. 
The vacuum terms indicate that the down-conversion 
process does not always produce photon pairs. In fact, 
down-conversion fails the vast majority of the time. 

If conditioning upon detection of a photon in one of 
the modes is performed, then upon success of the con- 
ditioning process there is high probability that a single 
photon is present in the other mode. The wave-packets of 
photons produced through conditioned down-conversion 
are characterized by a double-sided Lorentzian wave- 
function, 



2xk 



2ttx 2 k 2 



(21) 



where it has been assumed that the intrinsic response 
time of the conditioning detector Tdet, obeys 1/rdet 3> K. 
The frequency distribution of down-conversion sources 
has been extensively studied 0, EE EE EE 03 • 

By applying filtering to the conditioned mode, one can 
perform non-local pulse shaping [18L [l9j . An example of 
such a scheme is shown in FigurelH It has been shown by 
Ou E3 an( A Aichele et al. Eg that to obtain high purity 



of the post-selected state, narrowband filtering must be 
applied to the conditioned mode. Non-local pulse shap- 
ing is not limited to the temporal domain. In principle, 
shaping can be performed in any degree of freedom in 
which the photon pair are entangled. 

On the one hand, non-degenerate parametric down- 
conversion is quite suitable for quantum information pro- 
cessing applications due to its 'raw' symmetric profile and 
the ability for non-local pulse shaping techniques to al- 
low for approximately Gaussian (i.e. optimal) photons 
to be engineered. These properties have made spon- 
taneous down-conversion the system of choice for in- 
principle demonstrations in which efficiency is not an is- 
sue. However, although photons are heralded, scalability 
would require good quantum optical memories. While 
in- principle demonstrations have been performed |20j , ef- 
ficiencies are currently too low to be practical. 



C. Cavity QED pump pulse manipulation 

Perhaps the best solution is to combine a coherent 
excitation with single emitter technology. This can be 
achieved by using a Raman process to pump a single 
emitter in a high-Q cavity. 

Keller et al. [21j demonstrated that, through manip- 
ulation of the pump pulse, the temporal wave-function 
of photons emitted from Raman pumped single ions 
trapped in a cavity can be readily manipulated. Ex- 
perimental results using a single trapped 40 Ca + ion, 
demonstrated the production of Gaussian, rectangular 
and double-peaked pulse formations. 



V. CONCLUSION 

We considered the influence of photon wave-packet 
shape on the effects of mode-mismatch in linear opti- 
cal circuits, from which we established two criteria which 
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optimize circuit tolerance against such effects. Firstly, 
photons should be as broad as possible in the degrees of 
freedom in which mode-mismatch is likely to be intro- 
duced. Secondly, for a given bandwidth, photons with 
Gaussian profile are optimal. Our findings are completely 
general and hold for arbitrary linear optics circuits (i.e. 
beamsplitter networks). This includes ones which incor- 
porate post-selection and classical feed-forward, making 
the findings applicable to linear optics quantum comput- 
ing circuits. 

We considered various photon production techniques, 
discussing their advantages and limitations in producing 



photons suitable for quantum information processing ap- 
plications. Of these, some are inherently more suitable 
than others, within the context of our established crite- 
ria. 
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